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Momentum-space approach to calculation of one-electron energies and wave functions proposed initially 
by Fock for a hydrogen atom and considered later by Shibuya, Wulfman, and Koga for diatomic molecules 
is applied to clusters composed of three and more atoms. The corresponding basis set in the coordinate 
space is of the Sturmian type since all the hydrogenlike orbitals in this set have a common exponent, i.e., 
correspond to the same energy (as opposed to one-electron atomic orbitals). By the examples of He^ 7 and 
He^J" 11 cluster ions it is shown that increase in the number of orbitals in the set results in rapid convergence 
of eigenenergies and eigenfunctions of highly excited states. The momentum-space approach to the one- 
electron many-center problem may be used for various solid-state and quantum-chemical applications. 
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!>■ i I. INTRODUCTION 

<N ■ ....... 

In general, the one-electron many-center problem (i.e., the solution of the Schrodinger equation for an electron 

moving in the field of an arbitrary number of three-dimensional Coulomb centers) is known to have no exact analytical 

solutions. Making use of the variational approach does not allow for calculation of characteristics of excited states 

with controlled accuracy. Meanwhile, the knowledge of one-particle energy spectra and wave functions could facilitate 

the solution of numerous solid-state and quantum-chemical problems. 

£T) ' In the works fiLQ>H>Q> the nonvariational approach based on the Fock representation [j| of the momentum-space 

Schrodinger equation has been proposed as an alternative to the usual position-space methods. However, in the papers 

0, 0, EJ H the momentum-space approach has been applied to the diatomic molecules only (two-center problem) . 

The purpose of the present paper is to check the applicability of the momentum-space approach to the one-electron 

' many-center problem, i.e., to calculation of one-electron energies and wave functions (including those for highly excited 

states) of systems composed of three or more Coulomb centers. Among other things this is important for calculations 

on atomic clusters. It is worthwhile to note that the many-center problem is basically different from the two-center 

one in that there appear three-center overlap integrals which are absent in diatomics., and it is not obvious in advance 

' if the momentum-space approach will remain powerful in this case. 

Below we present the results of calculations of the one-electron energy spectra and wave functions of Hc^ 7 and 

He^" 11 cluster ions within the momentum-space approach. The results obtained are verified by the direct numerical 

solution of the coordinate-space Schrodinger equation on a grid. We use atomic units e = m = Ti = 1 throughout the 

paper. 
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II. MATHEMATICAL FRAMEWORK 



In a momentum space, the Schrodinger equation [—A/2 — E + V(r)]$(r) = for the energy E and the wave 
function 5'(r) of an electron moving in an external potential V(r) is an integral equation for a Fourier transform 
*(p) = (2tt)- 3 / 2 / dr*(r)exp(-ipr) of *(r) [1 0, II & : 

(|p| 2 + | Po | 2 )v|/(p) = -2(2^)- 3 / 2 J d V 'V{ V - p')tt(p') , (1) 

where |po| 2 = 2|_B| and V(p) = (2-7r)~ 3 / 2 J drV(r) exp(— ipr). For a cluster composed of Ni on ions with charges Zk 
and coordinates R&, where k = 1, Ni on is the number of an ion in the cluster, one has V(r) = — J^k -^Wl r — R-fc | 
so that V(p) = -47r(27r) -3 / 2 |p| -2 J2k exp(-ipR fe ), and hence Eq. Q reads 



|p| 2 + |p | 2 )*(p) = ^- 2 ]>> fc expHpR fe ) J dp'|p- P r 2 exp(zp'R fc )*(p') 



(2) 



The Fock transformation 5] projects the three-dimensional momentum space onto the four-dimensional sphere with 
a radius |po|, where |po| is the energy-related momentum that enters Eqs. QJ and The function ^(p) is related 
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to its four-dimensional image ^>(il) through 

*(p) = 4|p | 5 / 2 (|p| 2 + |p | 2 )- 2 *(f7) , (3) 

where Q stands for the set of three angles (a, 9, <f>) resulting from the Fock transformation of a momentum vector p = 
(|p|, 9, 4>). Here the angles 9 and <fi have usual meaning and the angle a is related to |p| through a = 2 arctan(|p|/|po|) 
so that < a < n for < p < oo. The function ^(fi) can be represented as 

tf(fi) = a kN exp(-ipR fc )y^(n) , (4) 

k,N 

where N — (n, I, m) stands for a set of three quantum numbers (principal n = 1, 2, orbital I = 0, 1, n — 1, and 
azimuthal m = —I, I) and Yat(JI) is a four-dimensional spherical harmonic defined as Yn(CI) = (—l) l C n i(a)Yi m (9, 4>). 
Here C n i are Gegenbauer polynomials and Yi m are usual three-dimensional spherical functions. 

By the example of a diatomic molecule, it was shown in Ref. 0] that coefficients OkN in Eq. 0} satisfy the set 
of coupled homogeneous algebraic equations. Generalizing the consideration of Ref. |3j to the case of an arbitrary 
number Ni on of ions in the cluster, we have 

^2H ijaj =0, (5) 

3 

where i = (k,N), j = (k',N'), 

Hij = |po|5^'(R fc - Rfe') - Zk " ^7^"( R fc - Rfc")S#»(Rfc" - Rfe') , (6) 

k" N" 71 

S%'(R k -%) = J dflY^(n)Y N ,{Q)exp[ip(K k - R fc ,)] . (7) 
The function S$ (Rfe — Rfc' ) can be expressed 000 as an integral 



S#'(R fc - H*0 = T7~T / rfr ^" (8) 
iPol J F — K]fe| 

over localized hydrogenlike orbitals 

XJV (r) = (-lr-'-^pol'/'Kn-i-lJIMn + l)!]^^ , (9) 

where i = |po||r| and L\ are the associated Laguerre polynomials. Note that all orbitals xn(t) have a common 
exponent |po|, contrary to one-electron atomic wave functions. The functions xn{y) form the Sturmian basis set 
which, as noticed in Ref. is a basis of a Sobolev space rather than a Hilbert space. All orbitals xjv(r) correspond 
to the same energy, regardless of their quantum numbers. 

Eq. JSJl is the nonvariational equation. To find the expansion coefficients a^N and the energy parameter |po|, one 
should solve the nonlinear equation det(Hij) = for |p |. Different solutions |p>o U °f this equation correspond to 
different eigenenergies Ei and eigenfunctions \&i(r) of the one-electron Schrodinger equation in the coordinate space. 
Once the values of |po| and akN are obtained, the corresponding electron energy and wave function are, respectively, 
E= -|p | 2 /2 and 

$(r) = ^a«w(r-Rt)- (10) 

k,N 

In general, there are four types of integrals over orbitals XJv( r ) m the matrix elements H^, see Eqs. (JSJ) and JSJ- 
First is a one-center integral 

dr X* N (r)XN<(r) (n) 
It appears if k = k', see Eq. (JSJ, and equals to (|po|/^)<5jvjv, so that S$ (0) = Snn 1 - Second is a two-center integral 
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X^r(r-Rfc)XiV'(r-R fe 
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that enters (Rfe — Rfc') at fc 7^ fc' (here R = R&/ — Rfe). Next, it can be shown 3] that at k — k' 7^ k" the sum 
53jv" m the second term of Eq. (JHJ equals to another two-center integral 




divided by |po|, where R = Rfe" — Rfe. Finally, at fc 7^ fc', fc 7^ fc", and fc' 7^ fc", the sum J2n" m ^1- © is a 
three-center integral 



divided by |p |, where Ri = Rfe< — Rfe and R 2 = Rfe" — Rfe. 

In the case of a single ion (N on = 1) with a charge Z, there are only one-center integrals (|llfl in the matrix 
elements H^, and Eq. © allows for a simple analytical solution [5j- Since iTjj = (|po| — Z/ti)5nn' sX k = k', one has 
|p | = Z/n and E = -Z 2 /2n 2 . 

As one can see from Eqs. © and ©, for a diatomic molecule (Nion = 2), the matrix elements Hij include both 
one-center and two-center Ijl2|l . (|13(l integrals. The analytical expressions for two-center integrals are rather 
complex 0. Hence, in order to solve Eq. JSJ numerically, one should restrict himself to a finite number M of orbitals 
Xn( t ) m the basis set. The authors of Ref. [3] applied the momentum-space approach to the molecular ion H^~. They 
calculated the energies of the ground (lsa g ) and one of the excited (2pa u ) states and showed that, as M increases, 
both those energies decrease and converge to the corresponding 'exact' values obtained by the analytical series solution 
in the coordinate space (see, e.g., Ref. |7|). 



If the number of ions in the system is greater than two (Ni on > 3), there appear three-center integrals 1)14(1 in the 
matrix elements Hij. In general, there are no exact analytical expressions for three-center integrals. However, because 
of the long-range nature of the Coulomb interaction, they are of the order of the corresponding two-center integrals. 
Hence, they cannot be neglected and have to be accounted for on an equal footing. 

In order to check the applicability of the momentum-space approach to the one-electron problem in the case Ni on > 
3, we have calculated the one-electron characteristics of several clusters with different ion charges and configurations. 
Below we present the results for chains composed of N ion = 4 and 6 ions with charge Z = 2 each. From the one- 
electron viewpoint, these systems correspond to cluster ions HeJ 7 and He^ 11 respectively. The linear form of the 
clusters was taken in order to facilitate the visualization of the wave functions in the coordinate space. We stress 
that our purpose was not to compute the energies and wave functions with an extremely high accuracy, as in Ref. Q 
for a diatomic molecule, but just (i) to clarify the very possibility to apply the momentum-space approach to atomic 
clusters, i.e., to the one-electron many-center problem and (ii) to study the convergence rate at the initial stage of 
increase in the number M of orbitals xw(r) in the basis set. So, we have restricted ourselves to M = 5, 14, and 23. 
These values of M correspond to account for all orbitals with n < 2; n < 3; and n < 4, 1 < 3 respectively. For each 
M, we have calculated the energies and the coordinate-space wave functions of 20 -j- 30 one-electron levels. 

In Table I, we list the one-electron energies for the first 15 levels of the chain He^" 7 . One can see that the energy of 
each level decreases with M, this decrease being more pronounced for highly excited states. As one goes from M = 5 
to M = 14, the energies of the lower four levels change by (0.3 -j- 0.5) %, while the energies of the higher levels change 
by (9 -i- 17) %. This is obviously due to larger weights of orbitals Xiv( r ) with n > 3 in the wave functions of highly 
excited states. Note, however, that further increase in M up to M = 23 results in an order of magnitude weaker 
change of the overall energy spectrum, less than by 0.015 % for the lower four levels and by (0.3 ~ 2) % for other 
levels. So, the convergence rate is very high, being comparable to that reported in Ref. [3J for H^" molecule. Roughly 
speaking, an increase in the principal quantum number n by one results in convergence of at least one more decimal 
digit, i.e., in at least an order of magnitude increase in the accuracy. 

By the direct numerical solution of the coordinate-space Schrodinger equation on a grid, we have verified that all 
eigenenergies considered indeed converge to their 'exact' values. This is true for the eigenfunctions as well. Figure 1 
shows the wave functions of the 5-th level calculated for different values of M . One can see that the wave functions 
for M = 5 and M = 14 differ considerably, while those for M = 14 and M — 23 are very close to each other. It is 
noteworthy that the wave function for M = 23 is visually indistinguishable from the 'exact' wave function obtained by 
the direct numerical integration of the Schrodinger equation. Note also that the small basis set (M = 5) is insufficient 
for even qualitative description of highly excited states. As one can see from Table I, in the case M = 5, the order 




(14) 



III. RESULTS AND DISCUSSION 



4 



of levels in the energy spectrum appears to be broken starting with the 14-th level. The 14-th and 15-th levels are 
degenerate for M = 5, the corresponding wave functions being equal to zero at the line connecting the ions in the 
chain. Meanwhile, those levels are non-degenerate for both M = 14 and M = 23, in accordance with the 'exact' 
numerical solution. Figure 2 shows the wave functions of the 14-th level calculated for M — 14 and M = 23. They 
differ considerably, while the wave function for M = 23 practically coincides with the 'exact' wave function. 

Figure 3 presents the one-electron energies of four lowest levels of the chain HeJ 7 as a function of the distance R23 
between two inner ions, the values of R12 and R34 being fixed at R12 = R34 = 2.4. The first and the second cluster 
levels emerge from the first levels of the diatomic He J 3 . As R23 decreases, the energy separation E 2 — E\ increases 
due to stronger hybridization between molecular orbitals. The same is true for the splitting E4 — E 3 between the 
fourth and the third cluster levels which both emerge from the second levels of the diatomic Hc^ 3 as they approach 
each other. Note that E 2 — E\ at a given i? 2 3 is much greater than E4 — E3. This is because of the different symmetry 
of the first and second diatomic orbitals (symmetrical and antisymmetrical, respectively). We have also solved the 
Schrodinger equation numerically at several values of R23 , making use of the finite difference method on a fine grid, 
and verified the results obtained by the momentum-space approach. 

Finally, we calculated the one-electron energies and wave functions for the first 29 levels of the chain Heg" 11 . The 
results are similar to those presented above for the chain He^ 7 . Increase in M leads to the progressive decrease of 
all eigenenergies and their rapid convergence to the corresponding 'exact' values. Figure 4 shows the wave functions 
of the 10-th level computed for different values of M. Again, as in the case of the chain He^ 7 , the wave functions 
for M = 5 and M = 14 differ considerably, while those for M = 14 and M = 23 are very close to each other, the 
wave function for M = 23 being almost identical with the 'exact' wave function. The eigenenergies are E10 — -2.414, 
-2.828, and -2.872 for M = 5, 14, and 23 respectively. 



IV. CONCLUSIONS 

In summary, we made use of the nonvariational momentum-space approach to calculate the one-electron energy 
spectra and wave functions of the ground and a large number of excited states of small atomic clusters. We have 
found that the one-electron characteristics converge rapidly with increase in the number of hydrogenlike orbitals in 
the basis set, each orbital having the same exponent, i.e., corresponding to the same energy. Our results show that 
the momentum-space approach to the one-electron many-center problem is rather powerful and may be considered as 
an interesting alternative to the position-space methods. 
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Table I. The one-electron energy spectra of the chain HeJ calculated for different numbers M of hydrogenlikc 
orbitals XJv(r) in the basis set. The value M = 5, 14, and 23 corresponds to the basis set composed of (1S,2S,2P), 
(1S,2S,2P,3S,3P,3D), and (1S,2S,2P,3S,3P,3D,4S,4P,4D) orbitals respectively. The distance between two inner ions is 
-R23 = 3, the distance between each outer ion and its neighbour is R12 — R34 — 2.4. The energies and lengths are 
measured in atomic units. 



Level number 


J\/f C 


A/f 1/1 


lvl — zo 


1 


-3.8961 


-3.9151 


-3.9157 


2 


-3.8460 


-3.8584 


-3.8586 


3 


-3.4430 


-3.4543 


-3.4543 


4 


-3.4348 


-3.4433 


-3.4434 


5 


-2.0881 


-2 3063 


-2.3149 


6 


-1.9496 


-2.2043 


-2.2237 


7 


-1.9496 


-2.2043 


-2.2237 


8 


-1.8304 


-2.0570 


-2.0678 


9 


-1.8304 


-2.0570 


-2.0678 


10 


-1.7869 


-1.9504 


-1.9848 


11 


-1.6536 


-1.9431 


-1.9532 


12 


-1.4881 


-1.7440 


-1.7712 


13 


-1.4881 


-1.7440 


-1.7712 


14 


-1.4177 


-1.6308 


-1.6657 


15 


-1.4177 


-1.5386 


-1.5621 
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Figure captions 

Fig. 1. The wave functions ^(x) of the 5-th one-electron level of the chain He^ 7 calculated for different numbers 
M of hydrogenlikc orbitals Xjv(r) in the basis set, M — 5 (dashed line), M = 14 (dotted line), and M = 23 (solid 
line). The ion coordinates are (-3.9, 0, 0), (-1.5, 0, 0), (1.5, 0, 0), and (3.9, 0, 0) in atomic units. For M = 23, the 
wave function practically coincides with that obtained by the direct numerical integration of the coordinate-space 
Schrodinger equation. 

Fig. 2. The same as in Fig.l, for the 14-th one-electron level of the chain He^ 7 . For M = 5, the wave function is 
zero at y = and z = 0, see the text. For M = 23, the wave function practically coincides with that obtained by the 
direct numerical integration of the coordinate-space Schrodinger equation. 

Fig. 3. The one-electron energies of the first four levels of the chain HeJ 7 versus the distance R23 between two inner 
ions, the values of R12 and R34 being fixed at 2.4 each. Solid lines are the results of the momentum-space approach 
for M = 23 orbitals in the basis set. Circles are the numerical solutions of the Schrodinger equation on a grid. 

Fig. 4. The wave functions \&(x) of the 10-th one-electron level of the chain He^ 11 calculated for different numbers 
M of hydrogenlikc orbitals XAf( r ) m the basis set, M = 5 (dashed line), M = 14 (dotted line), and M — 23 (solid line). 
The ion coordinates are (-6, 0, 0), (-4, 0, 0), (-1, 0, 0), (1, 0, 0), (4, 0, 0), and (6, 0, 0) in atomic units. For M = 23, 
the wave function practically coincides with that obtained by the direct numerical integration of the coordinate-space 
Schrodinger equation. 
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